Optomechanics with molecules in a strongly pumped ring cavity 
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Cavity coohng of an atom works best on a cychc optical transition in the strong coupling regime near res- 
onance, where small cavity photon numbers suffice for trapping and cooling. Due to the absence of closed 
transitions a straightforward application to molecules fails: optical pumping can lead the particle into uncou- 
pled states. An alternative operation in the far off-resonant regime generates only very slow cooling due to the 
reduced field-molecule coupling. We predict to overcome this by using a strongly driven ring-cavity operated in 
the sideband cooling regime. As in the optomechanical setups one takes advantage of a collectively enhanced 
field-molecule coupling strength using a large photon number. A linearized analytical treatment confirmed by 
full numerical quantum simulations predicts fast cooling despite the off-resonant small single molecule - single 
photon coupling. Even ground state cooling can be obtained by tuning the cavity field close to the Anti-stokes 
sideband for sufficiently high trapping frequency. Numerical simulations show quantum jumps of the molecules 
between the lowest two trapping levels, which can be be directly and continuously monitored via scattered light 
intensity detection. 



I. INTRODUCTION 

In recent years cavity cooling and manipulation of small 
polarizable particles has developed into an extended field of 
research 1 1, 2]. Considerable theoretical progress understand- 
ing the strongly coupled atom-field dynamics including par- 
ticle motion has been made in the past decades and 
has been accompanied by impressive experimental achieve- 
ments izl-lil. 

As recently discussed in a detailed review article ||1|], a gen- 
eralization to molecule cooling is limited by the inability of 
achieving sufficiently strong dispersive coupling to the cavity 
field. At the usable frequencies far from any internal molecu- 
lar resonance the ratio of polarizability and mass is orders of 
magnitude below the value close to resonance. Although cool- 
ing persists in principle it slows down too much to be useful. 
One way out of this dilemma is the use of collective enhance- 
ment using a high density of particles inside the mode vol- 
ume 10] . Again in principle the mechanism has been shown to 
work for pre-cooled atoms [13], but it turns out to be techni- 
cally challenging to implement the required high pressure and 
low temperature beam sources for molecules ^ . 

As an alternative route one can instead take a clue from the 
field of optomechanics, where motional degrees of freedom 
of large vibrating objects such as micro-mirrors lfT4|] . toroidal 
microresonators wall vibrations [15], nano-membranes |[T6h . 
levitated nanospheres [17] or even viruses Ell are manipu- 
lated via radiation pressure coupling to a light field. Such sys- 
tems not only exploit directed collective scattering (reflection) 
from a large number of very weak scatterers but in addition 
rely on generating a stronger light force just by applying very 
high light intensities. As the effective mechanical coupling 
grows with the intracavity field amplitude (i.e. the square root 
of the photon number), the inherently small single photon cou- 
pling rate, proportional to the extent of the objects zero-point- 
motion, can be augmented by several orders of magnitude to 
produce strong friction forces. 

In this work we consider the corresponding limit for a 




FIG. 1. A ring cavity with two degenerate modes (cosine: black line 
and sine: red line) with the same wave vector k. The cosine mode is 
strongly pumped and is responsible for trapping a polarizable particle 
around one of its intensity minima. Depicted in the figure is also 
the effective potential seen by the particle (blue line) as it is shifted 
around depending on the dynamical occupancy of the sine mode. 



weakly polarizable particle (molecule) in a ring cavity. In a 
previous analysis concentrating in the low photon number but 
strong atom-field coupling limit the ring cavity geometry has 
been shown to lead to enhanced cooling of atoms by only a 
factor two at best [ 19, 20]. Let us instead focus on the high- 
field limit for a symmetrically pumped ring cavity as in typ- 
ical optomechanical systems to study the corresponding cou- 
pling enhancement for a very weakly coupled molecule (see 
Fig.[T]). For symmetric pumping it proves useful to decompose 
the ring cavity field in a cosine and sine mode instead of the 
more common propagating mode expansion. In this basis pho- 
tons from the highly excited driven cosine mode are scattered 
into the almost empty sine mode. Amplitude and phase of 
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this weak single molecule light scattering to the sine mode is 
highly spatially dependent and strongly enhanced by the pres- 
ence of a large field amplitude in the cosine mode. The situa- 
tion for a high field seeking particle is illustrated in Fig.[T] The 
particle is trapped close to the antinodes of the highly popu- 
lated cosine mode, while the 7r/2-shifted sine mode provides 
the weakly-populated quantum mode that produces sideband 
cooling. Actually an analogous situation will occur for any 
cavity geometry with two almost degenerate modes with dif- 
fering mode functions. 

In the following using different theoretical approaches we 
will demonstrate two central physical effects for such a setup: 
i) the efficiency of cooling in the dynamical regime is consid- 
erably improved over cooling in a typical standing wave cavity 
with large circulating power and ii) for particles trapped in a 
standing wave cosine mode in a ring cavity the coupling of 
the single molecule to the corresponding sine mode is equiv- 
alent to an optomechanical coupling with a coupling strength 
enhanced by a power of the intracavity cosine mode photon 
number. 

The paper is structured as follows: in Sec. II we present an 
effective quantum Langevin description for the coupled dy- 
namics of a polarizable particle moving in the field of a ring 
cavity. In Sec. Ill we ignore noise terms and study impor- 
tant dynamical features of the underlying classical dynamics 
and compare this to the typical standing wave cavity cooling. 
In Sec. IV we concentrate on the case of an intense pump 
field with strongly confined particle motion, where the domi- 
nant interaction couples the sine mode to the particle motion. 
In Sec. V we numerically study the cooling process using 
wave function simulations and we present numerical exam- 
ples for quantum trajectories exhibiting the evolution of the 
molecule's state under continuous monitoring of the scattered 
field. Sec. VI contains more detailed discussions on the physi- 
cal interpretation of the cooling mechanism in this setup. Con- 
clusions are presented in the closing Sec. VII. 



II. QUANTUM LANGEVIN APPROACH TO RING CAVITY 
COOLING 

We consider a one-dimensional configuration as illustrated 
in Fig. [T] where a polarizable particle of mass m modeled as 
a two-level system with energy splitting hua is off-resonantly 
coupled to the field in a ring cavity. For a wave vector k (cor- 
responding frequency ujc) there are two degenerate counter- 
propagating modes of the cavity driven with amplitudes 7^+ = 
?7_ =77. One can rethink this in terms of cosine and sine 
standing wave modes, with the cosine mode driven with rj 
while the sine mode is unpumped. The two optical modes are 
represented by field operators ac (cosine mode), with commu- 
tator [ac,aj] = 1 and a (sine mode) with [a,a^] = 1. The 
cosine mode is driven by a laser at frequency ujp and ampli- 
tude T] with detuning A = ujp—ujc- Both fields decay at rate 2k, 
through the mirrors. The particle's quantized position and mo- 
mentum operators are x andp, [x^p] = ih. The field, strongly 
detuned from the particle's resonance by = 0;^ — ^p, is 
coupled with strength g so that we have very low saturation 



and the spontaneous emission 7 plays only a minor role, as 
discussed in more detail later in Sec.|Vll The position depen- 
dent intracavity intensity seen by a particle at position x can 
be written as (a J cos kx + al sin kx){ac cos kx + sin kx). 
Elimination of the internal degrees of freedom of the particle 
in the limit of very large detuning A^ ^ 7 leads to an effec- 
tive position modulated energy shift IJl 

U{x) = alaJJc {x) + alasUs{x) + (aja^ + acal) Ucs{x), 

where Uc {x) = Uq cos^ kx. Us {x) = Uq sin^ kxJJcsix) = 
Uq sin kx cos kx and 0^ g'^/ Aa. We arrive at an effective 
master equation for the system: 



where the Hamiltonian evolution is governed by 



(1) 



p 

H = hA [alac + alas) — hU {x) + ih (rjal — r]*ac) , 

(2) 

and the Liouvilian describes dissipation via cavity decay 
LcP = nD [ac] p and Lp = nD [a] p, where the super- 
operator's D action on a generic operator c is defined as 
D [c] = 2cpc^ — c^cp — pc^c. Notice that we have ignored 
spontaneous emission in the Liouvillian which amounts to two 
effects: decay of the excited level to the vacuum modes and 
subsequent momentum diffusion of the particle owing to ran- 
dom kicks produced by spontaneously emitted photons. The 
validity of our treatment will be analyzed in more detail in 
Sees. and IVll while here we make the observation that the 
scaling of the rates of these detrimental processes with A~^ 
makes them completely negligible for the typically huge de- 
tunings that molecules are driven at. 

We will in the following use the quantum Langevin equa- 
tions approach 



dacjt) 

dt 
dasjt) 

dt 

dx(t) p 

dt m' 
dp{t) _ _^dU{x) 



= [-K + iAc (x)] ac + iasUcs (x) + 7^ + «inX3a) 
= [-K + iAs {x)] as + iacUcs {x) + a-^, (3b) 

(3c) 
(3d) 



dt dx 

where the position dependent detunings are Ac,s (x) = A -\- 
Uc,s{x) and noise operators have vanishing correlation func 
tions except for ^ W 4 (O 

2kS {t-t'). 



III. CLASSICAL DYNAMICS 

Taking the quantum average of Eqs. ([3]) we get an infinite 
set of ordinary differential equations for expectation values of 
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FIG. 2. Cooling in standing wave cavity (blue) vs. ring cavity (red). 
Parameters [/o = 0.1 x At, Ac = —0.3 x n and 77 = 10.3 x 
(a) Field amplitude of cosine mode compared to the standing-wave 
cavity mode, (b) Sine mode field amplitude, (c) particle localization 
(d) momentum of particle. 



operator products. In a first approximation we can simply fac- 
torize any operator product expectation values into field and 
atomic operators and ignore noise operators, so that we ob- 
tain a finite closed set of corresponding equations for classical 
field modes and particles 1I21I] 

d , . 
d , . 



d , , 
d , , 



p_ 

m 



i<i + iAc {x)] ac + iasUcs {x) + r], (4a) 
+ iAs {x)] as + iacUcs {x) , (4b) 

(4c) 



= -h 



dU{x) 
dx 



(4d) 



The averages of cosine and sine modes field amplitudes are 
ac and as respectively, while the position and momentum ex- 
pectation value are x and p. First quantum correlations can be 
included in such a treatment by factorizing only higher order 
products as we will see later. 

Notice that the standing wave cavity model |3]can be re- 
produced by the above equations by setting Ucs (x) = 0. This 
enables us to perform a comparison between the efficiency 
of the ring cavity cooling and that of a standing wave cavity 
with equal driving 77. The numerical simulation results are 
presented in Fig. O The set of parameters that we consider 



= IMHz, A; = 27rxl0^ m" 



77a.m.u, Uo = 0.1 x a.:. 



T] = 10.3 X K, and initial momentum 5hk is unrealistic in that 
a high value of Uo around 100 kHz is feasible with atoms 
but not with molecules. The purpose of the exaggeration in 
the choice of Uq is to illustrate the advantage of ring-cavity 
cooling with large fields numerically since for smaller Uq in- 
tegration of Eqs. (O is a computationally hard and long task. 
The situation presented in Fig. [2] is that of a particle trapped 
inside a single potential well and cooled on a timescale of or- 



der 500 X K,~^ (as seen in Fig. |2]i) which is much faster than 
the corresponding standing- wave cavity cooling time. The de- 
tuning has been chosen at the value of A = —0.3 x k which 
is close to the trapping frequency of the particle in the cosine 
mode, choice that will be elucidated in the next Section. 



IV. THE STRONG CONFINEMENT LIMIT 

We will focus now on the quantum dynamics of this sys- 
tem in steady state for strong pump and deep trapping. We 
drop the index of ac and denote by a the steady state ampli- 
tude of the cosine mode, which we assume to be large and 
real (which can be guaranteed by the proper choice of phase 
of T]). For small co upling Uo <C we approximate it by 
a = \r]\ /V/^^ + A^. Note that as discussed in the following 
sections this only makes sense in the cooling regime, where 
we can expect a steady state state of the coupled system with 
a finite particle energy as well. 

In the limit a 1, one can ignore the term a^a sin^ kx in 
Eq. dJ]) as being small compared to the a'^ cos^ kx term, so 
that the cosine mode creates a deep periodic trapping poten- 
tial. Assuming good localization of the particle (kx <C 1) at 
the field maxima, one can expand cos^ kx 1 — k'^x'^ and 
sin 2kx 2kx. The Hamiltonian of Eq. ^ then becomes 



H 



P 

2m 



1 



2 ^2 
muj^x 



hAa^a - hUo{a + a^)x, (5) 



2 _ 



with effective harmonic trapping frequency a; 
2hUo{ka)^ /m and rescaled effective interaction strength 
Uq = Uoka. This interaction rescaling is similar to the case 
of optomechanical systems with mirrors or membranes in the 
limit of large photon numbers ll22l - [24ll where a fundamentally 
small coupling is collectively enhanced with the square root 
of the photon number. Note that a similar collective effect 
could be in principle used also for single standing wave mode 
cooling. However, as the atoms then are confined close to a 
maximum of this mode, where its derivative is zero we only 
get a coupling being second order in x, which is intrinsically 
much smaller leading to a weaker cooling effect. 



A. Fourier space analysis 

Reducing the operators to dimensionless quadratures, Q = 
x/xzpm and P = p/pzpm (where the ground state position 
uncertainty Xzpm = \/h/muJm and its corresponding mo- 
mentum spread Pzpm = ^/xzpm), one can write linearized 
Langevin equations in the form: 



— a(t) = {iA - K)a{t) + iUoQ + a^n, 
j^o){t) = {-iA - /^)at(t) - iUoQ + 4^ 

^P{t) = -oJmQ + Uoa + /7oa^ 



(6a) 
(6b) 
(6c) 
(6d) 
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The new coupling strength is 



(7) 



Note that the impHcit dependence of /7o on a through the 
Lamb-Dicke parameter kxzpm ~ ^m^^^ ~ reduces 
the collective enhancement scaling to Uq ~ a^/^. The differ- 
ence to the typical scaling with a as in optomechanical setups 
lies in the fact that the same field that mediates the interaction, 
namely the cosine mode, is also responsible with creating the 
optical trap. 

One can proceed to solve Eqs. Q as in 1I22I1 by deriving 
an effective reduced master equation for the motional degree 
of freedom or as in ||23|] by analytically deriving spectra of 
correlations of operators and integrating to find variances in 
steady state. For example, starting from Eqs. © one can 
derive the spectrum of the Q variance as {Q (uj) Q {uj')) = 
Sq (uj) S {cjO -\- uj') and integrate to find (AQ)^ n^t + 1/2. 
We state in the following results obtained via any of these 
methods in the perturbative limit where Uq <^ n. The physics 
of cooling in this regime is elucidated by considering the scat- 
tering rates into optical sidebands 



A{uj) 



An effective cooling rate is obtained as a difference between 
Antistokes and Stokes rates 



r A(uJm) - A{-UJm) 

and the final occupancy of 



nat 



/^^ + (cj^ + A)^ 
-4a;mA 



(8) 



(9) 



is optimized to ( / 2ujm ) under the condition of optimal cool- 
ing A = -Urn- 

The spectrum of the field fluctuations can be easily shown 
to be proportional to the spectrum of motion (cj) = 
A {uj) Sq {uj). Integration leads to a simple expression: 



A{uJm) A{-UJm) Up 

4A^r ~ -8A* 



(10) 



The apparent divergence in the expressions of n^t and is 
related to the fact that for very small detunings one leaves 
the cooling regime (l^ and the particle spreads out over a 
spatial range beyond the validity of the linearization. This 
behavior shows up as well in the numerical simulations of the 
full model (see Sec. |V]) and is more thoroughly analyzed in 
Sec. IIVBI where we show that the validity of our linearized 
treatment imposes that | A | cannot be smaller than the recoil 
frequency cj^ec = hk'^/2m. 



B. Validity of tight confinement treatment 



approximation imposes that kXzpm <C 1. In other words this 
requires that the trap frequency uom be much larger than the 
recoil frequency ujrec = hk'^ /2m, 

(jJrec ^ ^m- 

Moreover, localization of the particle in the final state de- 
scribed by the occupancy riat of Eq. Q requires that the 
spread of the final state (a thermal state) is also smaller than 
a wavelength, i.e, kxzpmy/'^at < 1- Assuming A small com- 
pared to (jJrn and ujrn of the order of k,, this condition is equiv- 
alent to 

|A| > U^rec 

As discussed in a number of papers ll22l - [24ll . sideband cool- 
ing works best when the Antistokes sideband is maximally en- 
hanced, i.e. when A = —uOrn- This implicitly sets a value for 



Moreover, to resolve the sidebands one has to also ask that 
ujrn > ^ which is equivalent to 



> 1. 



One can see that this requirement is not inconsistent with the 
assumption of our perturbative treatment, namely /7o <C a^, 
which, under the assumption that uJrn is of the order of k, can 
be reexpressed as 



Hi 



and which can be easily fulfilled for single molecules given 
the smallness of the coupling Uq- 
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FIG. 3. Atomic occupancy (a) and cooling rate (b) as a function 
of increasing 77 for fixed A = —2k. The minima of the full lines 
(i, ii and iii) occur at the point where the 3 dashed lines (that show 
— cJm/A) intersect the dashed black line of value 1, i.e. when the 
optimal sideband cooling condition is fulfilled. The cooling rates are 
also optimal, as expected, at values of 77 which ensure cjm = — A. 



The validity of the linearized treatment requires localization 
of the particle within an optical wavelength which in a first 
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C. Coupled equations for operator expectation values 

Starting from Eq. ([T]) where the Hnear dynamics is given 
by the simpHfied Hamiltonian of Eq. ([5]) and following the 
ideas shown in ll25ll one can derive a closed set of linear dif- 
ferential equations containing expectation values of pairwise 
operator products only. These can be solved dynamically and 
also exact steady state values for 

((at a), (Q2)^ (p2^) can be deduced: 



(a) 



(b) 



d 

~dt 
d 

di 



-2n{a^ a) - itjQ{{aQ) - (a^Q)), 



{Q^)=UJm{A), 



d_ 
di 

d 



dt 
d 

di 
d 

di 



(A) = 2uJm{{P^) - {Q')) + 

2Uo{{aQ) + (atQ)), 



(atp)), 



(11a) 
(lib) 

(11c) 
(lid) 



(aQ) = uJm{aP) - (a^ - iA){aQ) + iUo{Q^), (He) 
(aP) = {-K. + iA){aP) - uJm{aQ) + 



/7o ( (ata) + 1/2 + (a2) + -(^)), (llf) 



d 
di 



(llg) 



where {A) = {QP + PQ)- One can therefore derive steady- 
state solutions from this set of coupled equations for the mo- 
tional quadratures and intracavity photon number 



iP') = - 



{k^ + c^^ + A^){k^ + A^) + 2U^uJmA 
4u;^A(/^2 ^ + 8/72 A2 



(12) 
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(13) 

which in the limit where Uq <C are approximated by the 
expressions listed in Sec. lIV Al The method is a little more la- 
borious, but instructive since it provides one with results that 
hold for arbitrary coupling strength Uq and with some insight 
into the physics of the cooling process as seen from the rela- 
tion derived for the intracavity photon number 

This shows that is directly proportional to the imbalance 
between the uncertainties (Q^) and (P^), thus giving a good 
measure of the 'squashing ' of the uncertainty disk character- 
istic of a coherent state (or vacuum). 

Moreover, the exact expressions of (Q^) and (P^) allow us 
to make a clear investigation of the feasibility of our molecule 



FIG. 4. Log-log plots of atomic occupancy (a) and cooling rate (b) as 
a function of increasing 77 for variable A that is dynamically changed 
to fulfill the optimal cooling condition Um — —A. The occupancy 
goes to as expected for large 77 as the system goes deep into the 
resolved sideband regime (as seen in (a)). The cooling rate increases 
at first with 77^ and after entering the resolved sideband regime, its 
scaling drops to y/r]. The illustrations are for 3 values of Uq and 
clearly show that small final occupancies can still be obtained even 
with small Uq, however, at the expense of drastic increases in the 
cooling times. 



cooling scheme, i.e., into the cooling occupancy and rates 
with diminishing Uq. As seen in Fig.[3j for a fixed A = — 2^^, 
one can increase the cavity input power to reach a point of op- 
timal cooling where the occupancy is minimized as the cool- 
ing rate is maximized, and which is given by the condition 
00 rn = — A. This is valid for any arbitrarily small Uq. How- 
ever, as /7o is set to smaller values, the cooling rates are be- 
coming unrealistically small. Setting the optimal cooling con- 
dition, ujm = — A, the results of Fig. Ua) are straightfor- 
ward, showing that the occupancy asymptotically tends to 
with increasing r]. Investigating the analytical expression for 
r as given in Eq. ©, we see that in the unresolved sideband 
regime where ujrn ^ the scaling is with r]^ while deep 
in the resolved sideband regime where <C uJm the scaling 
drops to y^. This is clearly illustrated in Fig.^b) where the 
turning point occurs around an r] for which uJm ^ ^. 



V. NUMERICAL SIMULATIONS 

- After having demonstrated the key physical mechanisms 
present in the dynamics using different approximative ana- 
lytical treatments, we will now check the range of validity 
of these results by a numerical solution of the full quantum 
model. The quantum dynamics can be efficiently approxi- 
mated using time dependent Monte Carlo wavefunction sim- 
ulations of the full coupled atom field system. Here we can 
conveniently make use of our established previously devel- 
oped C++QED simulation package lH^ [27|] . Due to symme- 
try we can limit position space to one optical wavelength and 
truncate the photon numbers of the field modes sufficiently 
above the expected mean values, so that the overall size of the 
Hilbert space stays easily tractable. For very large pump pho- 
ton numbers, the cosine mode dynamics can be approximated 
by a coherent field, so that only two quantum degrees of free- 
dom have to be explicitly treated. This allows for even faster 
simulation and averaging over a large number of trajectories. 
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(a) (b) 




FIG. 5. For a), b) and c), parameters are Uo — 0.01 x n and ujm — 6k. a) Time dependence of average kinetic energy (mean over 100 
trajectories) of the particle in units of as a function of detuning, b) Cooling rate as a function of detuning, c) Simulated time averaged 
photon number in steady state as a function of detuning and coupling strength, d) Simulated and theoretical average sine mode photon number 
in steady state calculated in the linearized model as a function of detuning for Uq = k/20. 



A. Cooling rates and steady state temperature 

As a first step we confirm the efficiency of cooling as a 
function of cavity detuning by starting with a particle in a 
highly excited state and calculating its average energy loss 
as function of time for varying detuning between pump and 
cavity resonance. The kinetic energy as function of time is 
shown in Fig. Oa). Close to optimal detuning A of we see 
a fast decay towards the ground state kinetic energy hwm/^ 
with an optimal value of the detuning near A = —uJm- This 
value is different from the optimal value for /S. —k obtained 
for weak field cavity cooling and is characteristic for the op- 
tomechanical cooling limit where the mechanical oscillation 
frequency of the particle is bigger than the cavity linewidth 
K. Less cooling or even heating occurs for small detunings in 
contrast with the optomechanical model predictions; the rea- 
son is the breakdown of the tight confinement model as ex- 
plained in Sec. IIVBI A central important quantity for the 
usefulness of this cooling method is the corresponding cool- 
ing rate, which we plot in Fig. Ob). It also peaks near the 
sideband frequency uJm and reaches almost the order of mag- 
nitude of n well above /7o, which gives a practically quite use- 
ful timescale. 

Let us mention here that the simulations also give the steady 
state energy as a function of coupling strength and detun- 
ing. We find that for sufficiently large detuning we get ground 



state cooling also for very small coupling U^. However, as ex- 
pected for smaller Uq the timescale to reach the ground state 
gets longer. This is different from a standard optomechanics 
setup, where the mirror is additionally coupled to a thermal 
reservoir and the cooling rate itself largely determines the fi- 
nal temperature. The corresponding simulation results are in 
good agreement with the analytic results obtained above. 

Here we will show another aspect of this. As the cooling 
is connected to scattering of photons into the sine mode, its 
photon number directly provides information on the par- 
ticle's state and in particular on its temperature. The depen- 
dence of the sine-mode photon number on detuning and 
coupling strength is shown in Fig.Oc). It shows the predicted 
divergence at low detunings and agrees fairly well with the 
theoretical prediction from the linearized model above (See 
Fig. Od)), as in steady state the particle is well localized and 
linearization is justified. 

Let us note that that the sine-mode scattered photon num- 
ber, i.e. Ua is independent of mass, which is surprising ac- 
cording to a simple logic that a better localized particle would 
be expected to scatter less into the sine-mode, which is zero 
at its trapping point. We now add a direct comparison of the 
simulated and theoretically predicted average photon number 
in Fig. [3d). While for smaller detunings the agreement is ex- 
cellent, we see that at larger detunings the simulation predicts 
higher numbers. To a large extend this reflects the fact that 
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FIG. 6. Fixed parameters are Uo = /^/400 and A = —oJm- a) Averaged particle kinetic energy over 500 trajectories as a function of time 
for increasing trapping frequency starting with initial kinetic energy of Ek ~ 25hHi. b)Average steady state photon number in sine-mode. 
We see the decrease associated with the better localization in a deeper trapping potential, c) Dots give the cooling rate 7 averaged over 500 
trajectories. The solid line shows the expected power dependence derived from the linearized model above, d) Final steady state oscillator 
occupation number after a time nt = 300. 



even after a simulation time of uJrect = 1000 the system has 
not reached its steady state value. 

As next important issue we confirm the possibility of fast 
cooling even for very small coupling. For this purpose we 
show how the initial kinetic energy of a particle decays for 
small effective coupling Uo = a>:/400 with increasing pump 
power. A deeper potential leads to higher trapping frequen- 
cies, faster cooling and better final particle localization. This 
is visible in Fig. Oa), where we plot the loss of kinetic en- 
ergy of a trapped particle as an average over 500 trajectories. 
We see that for the same initial kinetic energy we get much 
faster relaxation with higher trapping frequency. This behav- 
ior gets even more obvious, if one looks at the average oc- 
cupation number of the trapped particle as function of time. 
While the final kinetic energy first drops due to better cooling 
but then increases with higher uorn due to the increasing zero 
point energy of the particle, the average particle vibrational 
excitation in the trap continuously drops and reaches almost 
zero for deep traps. Of course in practise, at some point one 
gets limited by the very small but yet nonzero spontaneous 
emission rate, which will induce heating and particle loss on 
a very long time scale. 

The cooling can be directly and nondestructively monitored 
in the scattered photon number shown Fig. Ob) which directly 
reflects better ground state cooling with increasing trap depth. 
The detailed correspondence of the two quantities as calcu- 
lated in the analytical section can be even more clearly seen 
in individual trajectories as discussed in the next section. 



One central question is the efficiency of the cooling for very 
far off resonance excitation, where the coupling of the parti- 
cle to the field modes gets small. Here we can make use of 
very high pump mode photon numbers as in optomechanics. 
Extracting the corresponding decay rate from an exponential 
fit in Fig.Oc) we see its nonlinear increase with 0;^ in accor- 
dance with the expected behavior form the linearized model. 

To show that this faster cooling does not enhance the diffu- 
sion significantly, we also plot the long time average of the os- 
cillator excitation number in Fig.Od), which predicts ground 
state cooling at large ujm- 

As a bottom line we clearly see that with sufficient power 
in the cosine mode to generate large mechanical trapping fre- 
quencies we can get fast molecular cooling despite very small 
coupling. As the induced photon number in the sine mode still 
stays small, we expect that this result holds even for more par- 
ticles cooled simultaneously as long the photon number in the 
sine mode does not increase too much. We plan to study this 
scaling in more detail in the future. 

B. Quantum Jumps 

One very useful advantage of the Monte Carlo simulation 
approach is that besides averages and expectation values one 
can also look at individual trajectories, from which one can 
get interesting clues on the underlying dynamics and direct 
hints for the expected experimental data. In Fig. [7] we show a 
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FIG. 7. Single trajectory expectation values showing the particle oc- 
cupation number, the photon number in the sine mode and the par- 
ticles position uncertainty as a function of time for A = —2k and 
other parameters as in Fig. [5] 

time window of the results of a typical trajectory as used above 
to find averages. In order to exhibit interesting dynamics we 
have chosen the case of a bit lower than optimum detuning 
A = —2k. < uom still yielding a non- vanishing occupation of 
the first excited trap level. In the simulation we continuously 
monitor the photons leaving the sine-mode. As we always 
know the momentary atom-field wavefunction, we can easily 
show its correlation with any other system variable. 

As expected we see a clear correspondence between the 
time evolution of the atomic position uncertainty and occupa- 
tion number and the average photon number in the sine mode. 
Note that after reaching the ground state the system shows 
striking and clear quantum jumps between the lowest and first 
excited state. While it spends long periods in the ground state 
we see sporadic jumps to the first excited state where the par- 
ticle remains for surprisingly long time until jumping back. 
This should be related with a symmetry change in the parti- 
cle wave function during such a jump. While such jumps are 
forbidden in a standing wave cavity setup, they attain a finite 
probability in a ring cavity |28], where spontaneous shifts of 
the trap minimum allow to break the mirror symmetry of the 
trap. As these fluctuations are rather rare, the first excited 
state, once excited, is much more long lived than the second 
excited trap state, which possesses the same symmetry as the 
ground state. The observation of such quantum jumps thus 
would clearly reveal the quantum nature of the molecule's 
center of mass dynamics. 

VI. DISCUSSION 

Lets now add some remarks on the effects of the spon- 
taneous photon scattering rate 7. Naturally a molecule is 
not a simple two-level system. Nevertheless considering a 
single ground state and one manifold of excited states cen- 



tered around uja with an energy spread bandwidth smaller 
than the laser detuning , can well describe the induced 
polarization dynamics. For weak saturation a two-level ap- 
proach with an effective dipole matrix element is then suffi- 
cient. In the limit of ^ 7, the induced molecular polar- 
ization can be decomposed into real and imaginary parts de- 
scribing dispersion (proportional IoUq) and absorption (pro- 
portional to 70 = 7^^/A^). To lowest order, 70 supplements 
the coupling strength Uq m Eqs. Q with an imaginary part 
^7o = ^ {kxzpm) X (7o<^)- Within our model, this implies a 
renormalizing of the coupling strength from /7o to \/^o+^- 
In view of 70 <C /7o, the effect is negligible. 

A second effect is the diffusion of the momentum triggered 
by the spontaneously emitted photons. A rate proportional 
to 7oQ^^ {kxzpm)'^ up to a geometrical factor is obtained that 
can again be neglected based on the smallness of 70 and the 
inhibition with the square of the Lamb-Dicke parameter. 

One can extend Eqs. Q to also include quantum fluctua- 
tions of the cosine mode. The fluctuations of this mode are 
close to those of the vacuum state with extra terms coming 
from the scattering of the sine mode photons off the molecule. 
Qualitatively this leads to a fluctuating trapping frequency oOm 
(that is proportional to a^) and a fluctuating coupling strength 
Uq. For a quantitative analysis, one can notice that in the 
regime which is optimal for cooling where A = —ujrn, the 
sine mode population is always close to zero (as shown in 
Eq. (fTOb ) and backscattering into the cosine mode can be ne- 
glected. When operating close to resonance, surely a large Uq 
can lead to a considerable number of photons in the sine mode 
and the back action and potential shift effects listed above 
have to be taken into account,^r example by following an 
approach similar to the one in ||3Q|] . 

As first noticed in (y\\\ , the linear optomechanical inter- 
action described by Eqs. Q can lead to entanglement be- 
tween the two degrees of freedom: vibration of molecule and 
quadratures of the sine mode, owing to quadrature-quadrature 
interaction. In terms of the logarithmic negativity, a maximum 
value of around 0.4 around the optimal cooling regime can be 
obtained and can be even surpassed by a proper filtering of 
the cavity output field as suggested in i29ll . As an advantage 
over the mirror-field entanglement scheme which is strongly 
affected by the thermal bath, the molecule-field entanglement 
is effectively obtained under zero temperature conditions. 

VII. CONCLUSIONS 

In conclusion we have shown that the ring cavity setup al- 
lows for strong coupling between a molecule's motion and 
the sine mode via scattering of sideband photons from the 
cosine mode into the sine mode. The low single molecule 
field coupling is enhanced by the large photon number of the 
strongly driven cosine mode. In this regime the dynamics is 
similar to the one described in typical optomechanical sys- 
tems with mirrors or membranes and we predict fast cooling 
rates. In contrast to superradiant cooling |0, [l3|] we do not 
need large molecular ensembles from start but rely on a single 
particle cooling mechanism. The central technical challenge 
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here seems to generate high enough intensities in a ring cav- 
ity to assure trap frequencies larger than the cavity Hnewidth 
even at very large pump detunings from any optical resonance. 
Here mirrors with a finesse over F = 1000 and supporting 
megawatt intracavity power are needed. From the molecular 
side one would like to have a high polarizability per mass with 
very low absorption at the chosen operating frequency. In this 
case even ground state cooling of any polarizable object seems 



feasible in a suitable ring cavity. As for atoms it could also be 
hoped that cooling properties get even more favorable in cav- 
ities invoking even more modes as in a confocal setup |2| . 
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